Coherence and Partial Coherence in Interacting Electron Systems 
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We study coherence of electron transport through interacting quantum dots and discuss the relation 
of the coherent part to the flux-sensitive conductance for three different types of Aharonov-Bohm 
interferometers. Contributions to transport in first and second order in the intrinsic line width of 
the dot levels are addressed in detail. We predict an asymmetry of the interference signal around 
resonance peaks as a consequence of incoherence associated with spin-flip processes. Furthermore, we 
show by strict calculation that first-order contributions can be partially or even fully coherent. This 
contrasts with the sequential-tunneling picture which describes first-order transport as a sequence 
of incoherent processes. 

PACS numbers: 73.23.Hk, 73.40. Gk, 73.20.Dx 



Introduction. The study of transport through quantum 
dots (QDs) revealed interesting phenomena like resonant 
tunneling, Coulomb blockade, and the Kondo effect. The 
measurement of the current, however, provides no infor- 
mation whether this transport occurs coherently or in- 
coherently. To approach this question the QD can be 
embedded in an Aharonov-Bohm (AB) geometry [Q-Q. 
A magnetic-flux sensitivity of the total current has been 
observed 0,^. Depending on the setup, the phase of 
the oscillations showed either a jump as a function of the 
gate voltage or a continuous phase shift. 

The transmission probability T^ ot (uj) through a QD for 
incoming electrons with energy oj and spin a is defined 
by its relation to the linear conductance, 
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In the absence of electron-electron interaction, trans- 
port can be described within a scattering approach |^,|^] 
with a transmission amplitude t^ ot (u)) oc GJ t c l R (w) and 



The Green's function G^ c lr(w) 
involves Fermi operators of the left and the right elec- 
tron reservoirs. It is, then, easy to show that transport 
through the QD is fully coherent. At low temperature it 
is possible to tune the transmission amplitude of a refer- 
ence arm and the AB flux such that the total transmission 
is zero, corresponding to a fully destructive interference. 
In the presence of electron-electron interaction, though, 
this approach fails. The transmission probability T^ ot (uj) 
can no longer be obtained from i^ ot (cj) introduced above 
but has to be determined using Green's function tech- 
niques for interacting systems Thus, the ques- 
tion of whether and how the coherent part of the trans- 
port through an interacting QD can be identified is non- 
trivial. It will be addressed in this Letter. 

First we use intuitive arguments to distinguish coher- 
ent from incoherent cotunneling through a noninteracting 
and an interacting QD. Then, for a quantitative analy- 



sis, we develop general expressions for the flux-sensitive 
transmission through an interferometer containing either 
one or two QDs. We derive the intuitively obvious result 
that the coherence of cotunneling may be spoiled by spin- 
flip processes. They give rise to an incoherence-induced 
asymmetry of the amplitude of the interference signal. 
We propose a symmetric AB interferometer using two 
QDs to show that first-order transport can be partially 
or even fully coherent, in contrast to the description of 
first-order transport within the language of incoherent 
sequential tunneling. 

We consider a single-level QD with level energy e, mea- 
sured from the Fermi energy of the leads. The Hamil- 
tonian H = Ht, + Hr + Hp + Ht contains H r = 
J2ka e kr%<rr a k<rr for tne left and right lead, r = LjR. 
The isolated dot is described by Hd = e J2 a n ° + ^ n T n i' 
where n a = c\c G , and H T = Y.kar^r^kar^ +H.c.) mod- 
els tunneling between dot and leads (we neglect the en- 
ergy dependence of the tunnel matrix elements tr,/ R ). 
Due to tunneling the dot level acquires a finite line 
width r = T L + T R with T L/R = 2n\t L/R \ 2 N L / R where 
Nl/r is the density of states in the leads. The electron- 
electron interaction is accounted for by the charging en- 
ergy U = 2Ec for double occupancy. To keep the discus- 
sion simple we choose U = for the noninteracting case 
and U = co for an interacting QD. 

Transport through noninteracting and interacting QDs. 
In the absence of interaction we can use the scattering 
formalism with the transmission amplitude 



(2) 



where the dot Green's function G^. et (o;) = l/(u> — e + 
iT/2) is the Fourier transform of -iQ(t)({c a (t), 4( )})- 
The transmission probability T^ ot {uj) = |^ ot (w)| 2 = 
FlTr/[(w — e) 2 + (r/2) 2 ] reflects resonant tunneling, 
which is fully coherent to all orders in T. 

We identify a contribution to the transport through a 
QD as coherent if by adding a reference trajectory fully 
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destructive interference can be achieved. Interaction of 
the dot electrons with an external bath (e.g. phonons) 
destroys coherence since interference with a reference tra- 
jectory is no longer possible: the transmitted electron has 
changed its state or, equivalently |L6]], a trace in the envi- 
ronment is left. Coherence can be also lost in interacting 
QDs by flipping both the spin of the transmitted electron 
and the QD. 

Away from resonance, |e| 3> ksT,T, transport is dom- 
inated by cotunneling fll3|,|Ii|]. There are three different 
types of cotunneling processes (for U = oo): 

(i) an electron enters the QD, leading to a virtual occu- 
pancy, and then leaves it to the other side. 

(ii) an electron leaves the QD, and an electron with the 
same spin enters. 

(iii) an electron leaves the QD, and an electron with op- 
posite spin enters. 

Process (iii) is elastic in the sense that the energy of the 
QD has not changed. It is incoherent, though, since the 
spin in the QD has flipped. 

The transmission, defined by Eq. (|l|) , through a single- 
level QD can be obtained ]l0|-|l2| from 
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For cotunneling, the transmission probabilities of elec- 
trons with energy ui near the Fermi level of the leads 
can also be obtained by calculating the transition rate in 
second-order perturbation theory and multiplying it with 
the probabilities P x to find the system in the correspond- 
ing initial state \. For an incoming electron with spin up 
the transmission probabilities are P X T lT /{Re[l / (uj — e + 
i0 + ) 2 ] with x = 0, T>i I0r case (i)> ( n )> an d (iii) , respec- 
tively. Since P + + p l = 1 and Po + P* = l/[l + /(e)] 
in equilibrium, where /(e) is the Fermi function, we find 
T$ ot {uj) = T^ ot ' coh (uj) +T^ incoh (Lu) with 



T^ ot (w) = Re 



r 7 T 
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We now show that Eq. (0) is not a good defini- 
tion for a transmission amplitude for interacting QDs. 
From 4 ot (uj) = i(P + Pf)y/TZT^/(u - e + i0+) we get 
\4 ot {uj)\ 2 = T^ ot (w)/[l + /(e)] 2 which not only does not 
yield the total transmission through the dot but also dif- 
fers from the coherent part of the transmission as well: 
there is no direct physical meaning of the expression 

|t*»| 2 . 

A generalization of the scattering approach has been 
proposed |18| ] which is compatible with the physical quan- 
tities expressed by Eqs. (|j) and (||). While this gener- 
alization is physically motivated, it gives no recipe how 
to calculate the transmission amplitudes explicitly in a 
given order in T. 



For [7 = 0, there are three more cotunneling processes. 
They involve double occupancy as an intermediate or ini- 
tial state. After summation, the spin-flip processes can- 
cel each other. In this case T^ ot ' coh (uj) = |4 ot M| 2 = 
T CT dot H =Re[^ L ^ /^ /(a,-e + ^0+) 2 ]. 

Inter •ferometry with a single QD. To support the results 
of our intuitive picture, we analyze quantitatively AB in- 
terferometers which contain one QD. The total transmis- 
sion probability T* ot (o;) through the AB interferometer is 
the sum of three parts: T^ ot {uj) and T rcf = |i rcf | 2 for the 
transmission through the dot and reference arm (the lat- 
ter is independent of energy u> and spin a), and the flux- 
dependent interference part T^ ux (ii;). Two kinds of ge- 
ometries have been considered, one using a two-terminal 
setup [j] and the other an open geometry . The two ge- 
ometries have in common that numerous channels (char- 
acterized by the energy u> and spin a) are probed simulta- 
neously, hence the interference signal is the sum of many 
contributions. To achieve fully destructive interference 
one needs to adjust the amplitude of the reference arm 
such that T^ ot (ui) = T lct for all contributing energies. 

We relate the flux-dependent linear conductance to the 
dot Green's function for the two-terminal geometry. To 
model the transmission through the reference arm we add 
to the Hamiltonian a term H IC { = J2kqa(t a kvR a q°'L + 
H.c.) with 2ttWN l N r = \t tef \e i( f. The AB flux $ enters 
via ip — 2ir$e/h (in a gauge that leaves the tunnel Hamil- 
tonian of the QD ^-independent). The current from the 
right lead is given by the time derivative of the electron 
number, / = ed(n R )/dt = i(e/h)([H,h R ]). The latter 
expression yields Green's functions which involve Fermi 
operators of the right lead. Using the Keldysh technique 
we relate these to the dot Green's function. After collect- 
ing all terms and using current conservation we find |l9| ] 
the surprisingly simple relation for linear response and 
first order in T and t (i.e., higher harmonics in ip are 
dropped) 



T flux, a(w) = 2v /r7r^ |i rcf | cos ip Re G^u;) • 



(6) 



For the second kind of interferometer it was shown ]?]] 
(under the condition that the open geometry ensures that 
the reference arm and the applied bias voltage do not 
affect the QD) that 



T« ux ' b H = 2 x /TLr^\f ei \Re[e- t0 GTH} 



(7) 



with 6 = ip + A9, where A8 is determined by the specifics 
of the interferometer. 

While Eqs. (§) and (|) are almost self-evident in the 
nonintcracting case, it was not a priori clear that they 
should hold for interacting systems as well. 

According to Eq. (^|) the conductance is always ex- 
tremal at ip = 0. Such a "phase locking" does not take 
place in the open-geometry setup: the AB phase at which 
the transmission is extremal can be continuously varied 
by tuning the energy of the dot level via a gate electrode. 
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In the absence of interaction the flux-sensitive interfer- 
ence part for the latter geometry is 



-ifiux.b . 



u) = 2\f ct \Re 
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At low temperature, ksT -C max{r, |e|}, we can adjust 
t Tei such that T^ ot (ui) — T rci for all contributing energies 
(up to corrections of order ksT/T and fcsT/|e|) which 
yields p| 



dP 



dV 



tx 4- 
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Here < 9 = arctan (r/2e) < ir. There is a value of 
the flux at which full destructive interference is achieved. 
When ksT > min{r, |e|}, the matching of all the trans- 
mission amplitudes does not work, and full destructive 
interference is not achieved. 

Let us now consider cotunncling (when |e| 3> T,ksT 
applies). Expansion of Eq. (0) leads to 
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showing that cotunncling in the noninteracting case is 
fully coherent. In the interacting case we find 
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dV 
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(11) 



For the two-terminal interferometer we obtain exactly the 
same as Eqs. ( |Io| ) and ( |lT| ) but with 9 replaced by <p and 
the proportionality replaced by the equality sign. 

The factor 1/[1 + /(e)] indicates an "interaction- 
induced" asymmetry associated with spin-flip cotunnel- 
ing, in accordance with Eq. (||). We, therefore, conclude 
that both kinds of interferometers discussed so far are 
suitable to distinguish coherent from incoherent cotun- 
neling through a QD. Moreover, this asymmetry is an ef- 
ficient and rather robust probe of the spin configuration 
of the QD (whether it has total spin or 1/2). The same 
information can be retrieved by the Kondo effect, how- 
ever, under more demanding experimental conditions. 

In most experimental situations, the number of dot 
levels participating in the transport exceeds 1. At a dis- 
tance |e| away from resonance, in the cotunneling regime, 
there are 2N relevant levels with iV~2|e|/A + l (here 
ksT, r <C |e| < Ec, the mean level spacing is A, and the 
factor 2 represents spin degeneracy). The number of lev- 
els within the range defined by temperature is 2M with 
M ~ 2fcsT/A + 1. The ratio of the number of coherent 
channels to the total number of transmission channels 
is (N + l)/(N + AM 2 ) in the valley where the electron 
number on the QD is even and N/(N + AM 2 ) when it 
is odd | p!9| ] . As a consequence the coherent contribution 
vanishes for fc^T 3> y/\e\A/8. Furthermore, the asym- 
metry between adjacent valleys diminishes for |e| ^> A/2. 



What about first-order transport, which dominates for 
fcfiT ^> r, |e|? The energy spread of electrons going 
through the reference arm is ksT, while the width of 
the resonance through the QD is T; hence, the matching 
of all the transmission amplitudes to the reference am- 
plitude does not work, and full destructive interference 
is not achieved. There is, however, at least partial coher- 
ence to lowest order in T. This manifestly contrasts with 
the sequential-tunneling picture which describes lowest- 
order transport as a sequence of incoherent tunneling pro- 
cesses. Thus, it does not take into account the coherence 
of the transmitted and reference beam, although it pro- 
duces the correct transmission probability through a QD 
in the absence of a reference arm. 

Interferometry with two QDs. The conceptual diffi- 
culty to address first-order transport in the above geome- 
tries is that the temperature has to be on the one hand 
large, yet, on the other hand, it has to be small to allow 
for a destructive interference of all energy components 
simultaneously. To circumvent this problem, we consider 
a two-terminal AB interferometer with two QDs, one in 
each arm. Then, fully destructive interference (in the ab- 
sence of interaction) is feasible at high temperatures. In 
related work, resonant tunneling (in the absence of in- 
teraction and flux) plj and cotunneling |^2|,^3| has been 
studied in the same geometry p3 . 

Each dot is described by the Hamiltonian introduced 
above for a single QD. We choose a completely symmetric 
geometry, and we assume ksT T, \ex\, | £2 1 as well as 
r ^S> |ei — ea|, where ei.2 is the energy of the level in QD 1 
and 2. In this regime lowest-order transport dominates, 
and we can set e = e% = £2- To model the enclosed 
flux we attach a phase factor e l(p ^ 4 to the tunnel matrix 
elements £_r,qdi and £_l,qd2, and e _lip / 4 to £l,qdi and 
£.R,QD2- The system is equivalent to one QD with two 
levels (each of them spin degenerate) with (^-dependent 
tunnel matrix elements. The total current is |10| 



with r L 



le 
2h 



dcjtr{[T L f L -T R f R ] G> 

[r L (i-/ i )-r*(i-/ 7? )]G<} (12) 

1 p +i<p/2 \ 



The matrices account for the two QDs. Expansion up to 
linear order in the transport voltage V and in the intrinsic 
line width T yields 



Qjtot 



dV 



Aire 2 



h 



du \f{u)Au{w) 



sin(yj/2) 



0G> 2 , M dG< 2 



+ [i-/H] 



(13) 



d(eV) Jy " d(eV) 

with Gf 2 M = Gf 2 (w) = 2iriP 2 1 5(oj - e) and Au(w) = 
S(u> -~ e) in the absence and An(u>) = S(cu — e)/[l + /(e)] 
in the presence of interaction. The off-diagonal density- 
matrix elements P 2 l = (|2)(1|^ vanish in equilibrium, 
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but they are present for finite bias voltages. To deter- 
mine them we use a real-time transport theory devel- 
oped in Ref. ]ll|,p^| and solve a generalized master equa- 
tion. We find [[L9| at V — and in zeroth order in T 
that 8P£/d(eV) = -(i/2)f'(e)shx((p/2) in the absence 
and dP£/d(eV) = -($/2)/'(e)/[l + /(e)] 3 sm(ip/2) in the 
presence of interaction. As a consequence, in the absence 
of an AB flux, only equilibrium Green's functions enter 
Eq. ([l3|). In the presence of flux, however, it is crucial 
to first account for finite voltage nonequilibrium Green's 
functions, and take the zero-bias limit only at the end. 
We find for the non-interacting case 



dP 



dV 



= 2- 



01 



dot 



V=0 



dV 



x [1 - sin 2 (^/2)] (14) 



with (dI dot /dV)\ v=0 = -(ire 2 /h)rf'(e) being the con- 
ductance through a single QD. At |sin(y>/2)| = 1 the 
total current vanishes, indicating that lowest-order trans- 
port is fully coherent. This completely contrasts the pic- 
ture of incoherent sequential tunneling. In the absence 
of interaction, however, the transport should be fully co- 
herent. For the simple limit U — we can rcdcrivc 
Eq. (|l4|) by using Eq. (7) of Ref. Q, determining the 
dot Green's function by an equation-of-motion approach, 
and expanding the result up to first order in T. 
In the presence of interaction we obtain 



OP 



dV 



dT 



dot 



v=o 



dV 



v=o 



sin 2 (^/2) 



[l + f(eW 



(15) 



with (dI dot /dV)\ v = = -(7re 2 // l )r/'(e)/[l + /(e)]. 

We point out that the total conductance is always 
smaller than the sum of the conductances through the 
QDs taken apart. The factor 1/[1 + /(e)] 2 yields an 
interaction-induced asymmetry in the ratio of coherent 
to total transport around a conductance peak. 

Conclusion. We have shown that interactions lead to 
an asymmetric suppression of destructive interference. In 
second-order transport we related this explicitly to spin- 
flip processes which give rise to an incoherent contribu- 
tion to the transmission probability. Even in first-order 
transport, the transmission is at least partially coherent. 
This statement is probably supported by the experiment 
of Yacoby et al. Q in which AB oscillations were observed 
in that regime. 

Our systematic analysis of how to describe the coher- 
ent components of physical observables in the presence 
of interaction (which is different from the way they can 
be accounted for in the absence of interaction) , may per- 
tain to other problems, such as the interpretation of the 
transmission phase through an AB interferometer [l], || . 
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